
#ifndef BTL_C_INTERFACE_BASE_H
#define BTL_C_INTERFACE_BASE_H

#include "utilities.h"
#include <vector>

template<class real>
class c_interface_base
{

  public:
	typedef real real_type;
	typedef std::vector<real> stl_vector;
	typedef std::vector<stl_vector> stl_matrix;

	typedef real* gene_matrix;
	typedef real* gene_vector;

	static void free_matrix(gene_matrix& A, int /*N*/) { delete[] A; }

	static void free_vector(gene_vector& B) { delete[] B; }

	static inline void matrix_from_stl(gene_matrix& A, stl_matrix& A_stl)
	{
		int N = A_stl.size();
		A = new real[N * N];
		for (int j = 0; j < N; j++)
			for (int i = 0; i < N; i++)
				A[i + N * j] = A_stl[j][i];
	}

	static inline void vector_from_stl(gene_vector& B, stl_vector& B_stl)
	{
		int N = B_stl.size();
		B = new real[N];
		for (int i = 0; i < N; i++)
			B[i] = B_stl[i];
	}

	static inline void vector_to_stl(gene_vector& B, stl_vector& B_stl)
	{
		int N = B_stl.size();
		for (int i = 0; i < N; i++)
			B_stl[i] = B[i];
	}

	static inline void matrix_to_stl(gene_matrix& A, stl_matrix& A_stl)
	{
		int N = A_stl.size();
		for (int j = 0; j < N; j++) {
			A_stl[j].resize(N);
			for (int i = 0; i < N; i++)
				A_stl[j][i] = A[i + N * j];
		}
	}

	static inline void copy_vector(const gene_vector& source, gene_vector& cible, int N)
	{
		for (int i = 0; i < N; i++)
			cible[i] = source[i];
	}

	static inline void copy_matrix(const gene_matrix& source, gene_matrix& cible, int N)
	{
		for (int j = 0; j < N; j++) {
			for (int i = 0; i < N; i++) {
				cible[i + N * j] = source[i + N * j];
			}
		}
	}
};

#endif
